The Potts model built on sand 
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Abstract: We consider the q = 4 Potts model on the square lattice with an ad- 
ditional hard-core nonlocal interaction. That interaction arises from the choice of 
the reference measure taken to be the uniform measure on the recurrent configu- 
rations for the abelian sandpile model. In that reference measure some correlation 
functions have a power-law decay. We investigate the low-temperature phase dia- 
gram and we prove the existence of a single stable phase with exponential decay 
of correlations. For all boundary conditions the density of 4 in the infinite volume 
limit goes to one as the temperature tends to zero. 



1. Model 

We define the model on the two-dimensional lattice Z 2 . Lattice sites 
are denoted by x, y, z and we write x ~ y if x and y are nearest neigh- 
bors. For a subset V C Z 2 we denote by dV the exterior boundary 
of V, i.e., the set of those y E Z 2 \ V such that there exists a nearest 
neighbor of y in V, V = VUdV, and the set d(Z 2 \V) is called the inner 
boundary of V. The square [— n, n] 2 fl Z 2 is denoted by V n . The con- 
figuration space is tt = {1, 2, 3, 4} z . Elements of are written as 77, £. 
For a configuration rj, i](x) G {1,2,3,4} is interpreted as the "number 
of sand grains" at x. For n 6 N, Q n = {1,2, 3,4} Vn denotes the set 
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of finite volume height configurations. Below we introduce measures 
H a p n on Q n that correspond to the finite volume Potts model at inverse 
temperature (3 with boundary condition a G {1, 2, 3, 4}, restricted to a 
special set of "recurrent configurations" defined from the abelian sand- 
pile model (cf. subsection 11.21 below). Our main result is that for (3 
large, n A B n forms the single stable phase of that model. 

1.1. Potts model. The Potts Hamiltonian with fixed boundary con- 
dition a G {1, 2, 3, 4} on the volume V n is 

x~y£V n 

That is a finite sum over nearest neighbor pairs of sites of which at 
least one belongs to V n and where it is understood that we substitute 
7](z) = a whenever z ^ V n . 

The Hamiltonians H n (rj\a) give rise to the finite volume Gibbs mea- 
sures on Q n : 

rf, n (v) = Jr exp[-/3# n (77|a)] 

where (3 > is the inverse temperature and where the normalizing 
factor Zg n is the partition function. It is well-known that there exists 
a critical inverse temperature (3 C G (0, +oo) such that for (3 < f3 c the 
Potts model has a unique infinite volume Gibbs measure (as n j +oo), 
i.e., the weak limits of v% n for different a coincide, while for (3 < (3 C , 
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these weak limits are all different and define 4 mutually singular ergodic 
Gibbs measures z/§ on fl, called the pure phases. 

1.2. Sandpile model. The abelian sandpile model in volume V n is a 
Markov chain on Q n . We briefly introduce that Markov chain, more 
details can be found in the original paper [[}. and in [2] or 0. 
The Markov chain starting from rj £ Q n is defined as follows. Suppose 
that rjt-i is the configuration at time t — 1 > 0. Pick randomly a site 
in V n , say x £ V n , and add one grain at x to r\ t -\. In case 7ft_i(:r) < 3, 
the new configuration is simply 

Vt(y) = Vt-i(y) +$x, y 

with S x>y the Kronecker delta. In case f]t-i{x) = 4, by adding one 
grain at x the number of grains at x becomes equal to 5. That site 
will now topple, i.e., 4 grains are removed from x and one grain is 
given to each neighbor of x in V n . At the boundary, grains are lost 
when the site topples. It is now possible that the number of grains 
at one or more neighbors of x exceeds 4 and we have to repeat the 
toppling operation on all of these, and so on. It turns out that no 
matter in what order we perform these toppling operations, at the end 
of the avalanche a unique configuration i] t £ £l n appears. In that way, a 
discrete time Markov chain on Q n is defined where the only randomness 
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is in the independently repeated uniform choice of the site where a grain 
is added. 

Analysis of that Markov chain learns that it has a unique class TZ n 
of recurrent configurations and the stationary measure A n is uniform 
on that class: 

I '*"n I 

see 0. 

Whether a particular configuration r] E Q n belongs to 1Z n can be 
decided from the output of the so-called burning algorithm |2]. The 
burning algorithm has as an input the configuration rj and its output 
is a set A C V n . It runs as follows: start from Aq = V n and remove 
("burn") all those vertices x G A (and edges containing x) which 
satisfy rj(x) > ua ( x ) where n v (x) denotes the number of neighbors of 
x in V. This gives A\\ now proceed in the same way with Ax, etc. until 
no further vertices can be removed. The output A of the algorithm 
is the set of remaining vertices. Recurrence is then characterized by 
"burnability" , i.e., r\ G TZ n if and only if A = 0, i.e., all vertices can be 
burned. 

The stationary measure A n is thus the uniform probability measure 
on all burnable configurations in V n . The cardinality \7Zv\ (= the 
number of recurrent configurations in V) equals the determinant of the 
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discrete Laplacian on V with open boundary conditions see . E.g. if 
V is a square containing N sites, then \TZy \ — (3.21)^. 

For a proof of these facts, see e.g. 0, jZjj or [T2]. Remark that for 
all finite V C Z 2 , the constant configurations r) = 4 and 77 = 3 are in 
IZy, but 77 = 2 and 77 = 1 are not recurrent except for some very special 
choices of V. One easily concludes that the condition that 77 G 1Z n is a 
nonlocal hard-core constraint. 

The following proposition is an immediate consequence of the burn- 
ing algorithm. 

Proposition 1.1. If 77 G TZ n and £ > r\ (pointwise), then £ G 7^ n . 

In some aspects the abelian sandpile measure X n , n — > 00, behaves 
as a model of statistical mechanics at the critical point, a phenomenon 
which is sometimes referred to as "self-organized criticality" because 
there is no explicit tuning of parameters. In the physics literature vari- 
ous critical exponents related to the avalanche behavior are introduced 
for that model. One signature of "critical behavior" is the presence of 
power law decay of correlations for the height 1 two-point function, as 
proven by Majumdar and Dhar in [I]: 

Proposition 1.2. There exist constants c,C > such that 

c\x\- 4 < \X n (v(0) = V(x) = 1) - A n (77(0) = l)X n (v(x) = 1)| < C |x|" 4 
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for all x 7^ and n large enough. 

On the other hand, a contour of 4's completely decouples the inside 
and the outside, as we now show. A subset V C Z 2 is called simply 
connected if the corresponding t^cK 2 obtained by "filling the squares" 
of V is simply connected. 

Proposition 1.3. For W C V n denote by A w the event that r)(x) = 4 
on W . For any simply connected subset V C Z 2 with V C V n : 

^n(VWv n \v\ 4 dv) = K(Vv\4 d v)K(Vv n \v\ 4 dv) (1-1) 

Proof. Denote by TZy\y the set of configurations which are burnable 
in V n \ V and such that the extension rjy^yAy is burnable in V n . By 
the burning algorithm, rjv4:dvVv n \v e ^v„ ^ an< l om y if Vv G T^v an d 
r] VA y E R-y^y Therefore 

I[riv e Kv,r) VnW e Uf 

text 

-v n \v\ 



MVWv n \v\ A 9v) = m^-WKvl 



KiwHdv) = r^H (1-3) 

I '^v\ 

and 

Uv VnW M = ^ ,^ y , (1-4) 

which gives the result. □ 
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1.3. Sandpile model with Potts interaction. Define the probabil- 
ity measures //§ n on fl n as 



a i \ exp(-(3H n (r]\a))I[r] e K T 



where the normalizing constant 5 ~ n is 



E ln = E exp(-/3iy n (77|a)) 



Similarly, we define the partition function v in an arbitrary finite 
volume V. This partition function will also be abbreviated as H 0j y. 
//^ n is of course just the original Potts measure conditioned on being 
recurrent: 

Obviously, at infinite temperature, (3 = 0, we recover the station- 
ary measure A n of the sandpile model. The constraint r\ 6 1Z n can 
be viewed as introducing an extra nonlocal hard-core interaction (im- 
plicitly given by the burning algorithm) but it also breaks the Potts- 
symmetry: approximately for n | +oo, \ n (r](0) = 4) = 0.4, X n (r](0) = 
3) = 0.3,A„(r/(0) = 2) = 0.2, A n (^(0) = 1) = 0.1 as computed by 
Priezzhev, |10j . 
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2. Results 

With boundary condition a = 4, at low temperature, the typical 
configurations of the Potts model on sand look like an ocean of 4's 
with exponentially damped burnable islands. 

Theorem 2.1. For any e > there exists (3q G (0, oo) such that for all 
(3 > P and all n G N 

A4,„fo(0) = 4) > 1 - e (2.1) 

Moreover there exists c > stzc/i £/iat /or (3 > (3q and n big enough we 
have the bound 

exponentially small in the distance \x\ from the origin. 

()2.2j) must be contrasted with the situation for f3 — where there 
are long range correlations, see Proposition 11.21 

Besides "all 4", the "all 3" is the only other groundstate. But that 
one is unstable: 

Theorem 2.2. For every a > there exists /3(a) G (0, oo) and c = 

c(a, (3) > such that for all [3 > (3(a) 

4M X e V n ■ V(x) = 3}| > a\V n \) < e- c l y "l (2.3) 
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for large n. 

The next section gives the proof of Theorem 12.11 and introduces a 
random cluster representation of the Potts model on sand. Section 0] 
is devoted to the proof of Theorem 12.21 It will be seen that, as an 
extension of Theorem I2.2[ for no matter what boundary conditions, 
the density of 4 tends to one with (3 "J* +oo. 

3. Random cluster representation 

The volume V n = V n U dV n can be considered as a finite graph with 
the sites x e V n U dV n as vertices and with edge set B n = B consisting 
of the nearest neighbor bonds x ~ y where at least one neighbor is in 
V n . We define the sand-Potts random cluster measure n = ip® on this 
graph with parameter p e [0, 1] as the probability measure on {0, 1} B 
which to each a e {0, 1} B assigns probability 

(pfa) = — [JJp <T(e) (l-p) 1_ff(e) ] Yl ^ is constant on clusters of a] 

P eeB n&Zn 

(3.1) 

By cluster we mean a (nearest neighbor) connected component of sites 
(including isolated sites) as obtained from the bond configuration a. 
Bonds for which a(e) = 1, a(e) = are called open, respectively, closed. 
In this definition, we assume that the boundary sites are all connected 
(wired). All sites that are connected to the boundary are in the same 
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cluster. The restriction that rj is constant on clusters also implies that 
i] is constant equal to a on the cluster of the boundary. Remember 
however that r\ = 2 and r] = 1 are not in 7£ n . 

3.1. Stochastic domination. 

Lemma 3.1. Lei a = 3 or a = 4. For every edge e = (xy) in B and 
every a B \ e G {0, l} B ^ e \ 

^(o-(e) = l\a BV ) = p (3.2) 

i/x and y are connected via open edges in o~B\e- If, on the other hand, 
x and y are not connected via open edges in o- B \ e , then we still have 

P>V a Me) = l\o- B \e) (3.3) 

and for a = 4, 

#) = 1W>^ (3.4) 

Proof. Let a G {0, 1} B . We write 

7[?7 is constant on clusters] = k(n, a; a) (3.5) 

for the number of recurrent configurations that are constant on the 
o~ — clusters and fixed equal to a for each site that is a— connected 
to dV n . It equals \R n \ when all edges in a are closed. Obviously, 
k(n, a; a) < \R n \ and k(n, a; a) is decreasing in a and is increasing in a. 
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For a = 3 or 4, k(n, a; a) > 1. Continuing with either a = 3 or a = 4, 
we have 

^p(l e CTB\ e ) pk(n,a; l e a B \ e ) 



<p%{O e a B \ e ) (1 - p) k(n, a, e a B \ e ) 

and hence 



(3.6) 



<$(a(e) = l\a BV ) = = ( 3 - 7 ) 

"T (pg(l««T B \ e ) ^"^V 1 ^ fc(n,a;l et T fl \ e ) 

Abbreviate cr 0,e = O e cx j B\ e and cr 1,e = l e crs\ e ; they are both equal to 
<7B\ e off e and <7°' e (e) = and a 1,e (e) = 1. 

To prove the first statement ()3.2|) : suppose x, y are connected via open 
edges in a B \ e , then every configuration r\ G lZ n compatible with a has 
r](x) = r](y), and hence k(n, a; a) does not depend on <r(e) in that case. 
However, if x and y are not connected via open edges in a B \ e , then 
we must investigate the effect of merging two clusters. By making e 
open, we connect two clusters and we must estimate the new number 
of recurrent configurations that are constant on clusters in terms of the 
old. Since always 

k(n, a; a°' e ) > k(n, a; a 1 ' 6 ) (3.8) 

we obtain (j3.3|) from ()3.7|) . For the last statement ()3.4|) . we combine 
Proposition 11.11 with (|3.7J) . Suppose that r\ 6 7£„ and is constant on 
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clusters C\ and C2 taking there the values a\ and a 2 respectively. The 
new configuration £ defined as 

£(x) = r)(x), x C\ U C 2 , £(2) = max{ai, a 2 }, x G C\ U C 2 

is still recurrent and is constant on Ci U C 2 . Moreover, if say C\ is the 
boundary cluster, then necessarily a\ — 4 and hence also max{ai, a 2 } = 
4 remains compatible with the boundary. (This does not work with the 
boundary condition a = 3.) Simple counting shows that the map r\ — > £ 
is at most seven to one, or 

fc(n,4;a°' e ) < 7Jfe(n,4, <r 1,e ) (3.9) 

Combination of O and flH2j) gives (Q- □ 

Let ^ g be the Bernoulli product measure on {0, 1} B with density 
q = ip q (a e = 1). 

Proposition 3.2. The random cluster measure ip 1 ^ stochastically dom- 
inates ipq with q = p/(7 — 6p), i.e., (Pp(o~(e) = l,e G E) > i[) q (a(e) = 
l,e G E) for all edge sets E. For a = 3,4, ipp 1 is stochastically dom- 
inated by ip p . Finally, (pp 1 always stochastically dominates (p a , for < 
p' < p < 1. 
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Proof. The first and second statement follows directly from Lemma l3~T| 
see e.g. Theorem 4.8 in j^j. The last statement follows from (|3.7jl and 



3.2. Coupling. The previous construction and arguments are analo- 
gous to and inspired by the Fortuin-Kasteleyn representation of the 
standard Potts model. To recover the q— state Potts model, one should 
simply replace 1Z n in ()3.5|) with Q n . Our next step, making a coupling 
between the 77— and the a— field, is the analogue of the Swendsen- 
Wang-Edwards-Sokal coupling, jTSlE]- For a general reference, see 0. 

We make a coupling ¥p n = between the Potts model on sand 
and the sand-Potts random cluster measure. Let P° be the probability 
measure on Q n x {0, 1} B constructed as follows. Assign first to each 
site in V n a sandvalue according to the probability measure A n and 
each site at the boundary dV n gets the value a. Independently, let 
each edge in B take the value or 1 with probabilities 1 — p and p 
respectively. Secondly, condition on the event that no two neighboring 
sites (including sites at the boundary) with different heights have an 
open edge connecting them. In a formula, 



(ED. 



□ 




e={xy)aB 
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where in the last indicator function it is understood that T](z) = a for 
z G dV n . 

Proposition 3.3. Suppose (3 = — ln(l— p). Then, [i% n is the marginal 
ofFp projected on fl n and ipp is the marginal ofF® projected on {0, 1} B . 

Proof. The proof is by direct computation. For example, if we sum 
over the a we have 

II [P u{e \ l - PY~ u{e) tMx) - r]{y))a{e) = 0]] = (1 - p yiv(*)?v(v)] 

a eG-B 

which determines 1 — p = exp[— (3\. □ 

3.3. Proof of Theorem I2.lt stability of the 4-phase. We apply 
Proposition 13.31 

= 4 ) 

> Pp(?7(0) = 4| is in the cluster of the boundary) x 

V?p(0 is in the cluster of the boundary) 
= (f p {0 is in the cluster of the boundary) (3.10) 

and now Proposition 13.21 to conclude that 

^ n (^(0) = 4) > V ; p/(7-6p)(0 is connected to the boundary) 
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which goes to one, uniformly in n as (3 = — ln(l — p) goes to +00. 
That shows (|2.1|) . For the exponential decay of correlations, ()2.2|) . 
observe that by the very same argument as above one shows that for 
large (3, in ^ n there is percolation of 4's uniformly in n. For each 
n the jJLp n — probability that there is a nearest-neighbor path of 4's 
connecting the origin with the boundary dV n is not smaller than the 
percolation probability in the Bernoulli bond process with occupation 
probability (1 — exp[— /3])/(l + 6 exp [—/?]). Moreover, we can always 
consider a rectangle parallel to V n between the origin and site x with 
one side proportional to \x\ and the other side equal to n. Again 
by the same domination argument, the /ii n — probability that in that 
rectangle, there is percolation of 4's from one side of V n to the opposite 
side is not smaller than 1 — exp[— c\x\] with c f +00 as f3 \ +00, 
uniformly in n. Therefore, denoting that event by "crossing', we have 
that 

l/^nfa^M )) -^A„(^(^)^( )l cr °ssing)| < exp[-c|x|] 

That can be combined with Proposition II if the origin and the site 
x are separated by a path of 4's, then they are independent: 



vi,n(v(x)v(Q) (crossing) = fj.% n (rj(x) |crossing)/ij n (^(0) (crossing) 
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For each of the two factors we can use the previous argument to con- 
clude the proof. 

4. Instability of the 3-phase 
In this section we prove Theorem 12.21 

The main idea is to consider the restricted ensemble defined below 
(following ideas from [3j. Let V be a finite (large) volume and r\ e fly 
a configuration in the volume V. Define the contours of 77 as connected 
components of edges separating the different values of 77 in neighboring 
sites (we now consider these sites as centers of unit squares and the 
edges are the sides of these squares). We suppose that the configuration 
77 G IZy is burnable. The set of configurations that have contours with 
lengths not exceeding 12 not touching dV is denoted by IZy . That set 
of configurations defines a restricted ensemble. It is easy to see that 
we have the following alternative: either the configuration rj is a sea of 
4 with islands made from contours with length < 12, or it IS cL SGcL of 3 
with such islands: 

7^ 2 = 7^ 2 4 U 7^ (4.1) 

Notice that of 2 or 1 is not possible by burnability. 

We list all sub configurations of energy < 12 in the sea of 4 and in the 
sea of 3. Of course, these islands have to be burnable. In the figure, 



empty cells are part of the sea. 
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4.1. Islands in the sea of 4. 



;i) Energy =4: 



(2) Energy =6: 



and rotations and reflections 



(this remark will be omitted from here on) 
(3) Energy = 7 



(4) Energy =8 

where a, b = 1, 2, 3. 
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b 



(5) Energy =9 



where a = 



a 


b 




b 



,2,3, b = 2,3, a ^ b. 



Energy =10 



a 


b 




c 



where a ^ b, b ^ c; a, 6, c = 1,2,3 and if a = c = 1 then 
6 = 3. 



where a = 2, 3, b = 1, 2, 3. 
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a 


a 


a 



a 




a 


a 




a 



where a = 2, 3 



Energy =11 



where a,b,c— 1, 2,3, a ^ 6. 



where b — 2, 3, a, c = 1, 2, 3, a ^ 6, a ^ c, 6 ^ c. 



a 


a 


a 




b 
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b 



where a = 2, 3, 6 = 1,2, 3, 6^ a. 
, where a=l,2. 
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(8) Energy =12 
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, where a,b,c = 1, 2, 3. 



where a, b, c — 1, 2, 3, a ^ b, a ^ c, if a — 1 then 6 = c = 3, if 
a = 2, then 6 = 3 or c = 3. 
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where a=2,3, b, c — 1, 2, 3, a 7^ b, a ^ c; if a — 2, then 6 = 3 
or c = 3. 
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for a = 1, 2, 3. 
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where a = 2, 3 and 6 = 1,2, 3. 
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where a = 2, 3. 



4.2. Islands in the sea of 3. 

Lemma 4.1. Islands with energy < 11 in the sea of 3 are in one- 
to-one correspondence with islands with the same energy of 4- The 
correspondence is given by substitution 3 < — > 4. 



Lemma 4.2. The same correspondence can be applied to all islands 
with energy =12, with two exceptions: 



1 




1 




1 





and 
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2 


2 




2 


2 



The correspondence gives all islands in the sea of 3 with energy 12. 



Proof. This follows from the fact that 
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and 



are burnable but 
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and 



are not burnable. 
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2 



□ 



i 712*) i ,(12) 

4.3. Free energy estimates. Define S 4 y and S3 y to be the partition 
functions of the ensembles 1Zl 2 v , respectively 1Z^ V . 



Lemma 4.3. There exists j3 such that for any f3 > /3 there exists a 
constant c > such that for any sufficiently large volume V 



logSi;?-logS^> C |T/|e-^ 



(4.2) 



Proof. Define the extended ensemble TZy^ xt , consisting of all configu- 
rations of disjoint burnable islands with energy < 12 in the sea of 3 in 
volume V. Evidently 



7? 12 r~ 7? 

and the inclusion is strict because e.g. 



12,ext 
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2 
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is a subconfiguration allowed in the ensemble TZy'^ xt but not in the 
ensemble 1Zy 3 . Indeed, it is not burnable but the island 



is burnable. 



If we analogously define TZ VA X , then in fact lZ v f x = 1Zy A because 
a configuration of (disjoint) islands in a sea of 4 is burnable if and only 
if all the individual islands are burnable in the sea of 4, see Proposition 
PI 

Therefore 

^(12,ext) . ,-,(12) 
"3,V — ^3,V 

and we can obtain Lemma 14.31 from the following 
Proposition 4.4. 

logH^-logH^^clyle- 1 ^ (4.3) 

The proof is an application of the usual polymer expansion (see 
(HI E]) because islands in the sea of 4 do not interact (again from 
Proposition ll.3|) . Comparing the expansions of log E^y and log E^y 
term by term, and using Lemma [4. II and Lemma \4.2\ we obtain a dif- 
ference in the terms of order e~ 12/3 and no difference in previous terms. 
In the same way it is easy to prove a weaker inequality valid for all 
finite V C Z d 
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Proposition 4.5. For (3 > (3 there exist c > 0, f G R such that for 
any finite V C Z d : 

log 2$ -logSgJ > logS$ -togSgJ^ > C |y|e-^-/i^| (4.4) 

□ 

4.4. Proof of Theorem 12. 2L Again, the connected components of 
the set of lattice edges separating different values of neighboring spins 
are called contours. Each contour 7 will be associated with values of 
spins on all squares touching 7. Consider the set F = . . . ,7„} of 
all those contours of a given configuration 77 e 7Zv,3 which have energy 
(length) exceeding 12, > 12. We call T the polycontour of r\. The 
polycontour separates the volume V into a finite number of connected 
domains V\, . . . , V m . To each domain Vi is associated the value of the 
spin Ki, induced by contours jj neighboring this domain (so K{ — k if 
the spins in the inner boundary of Vi are equal to k). If T = then 
there is the domain V, to which is associated k = 3. Denote the band 
consisting of squares touching T by [r]. From burnability it follows 
that all domains Vi for which Ki = 1, 2 are contained in [r]. We denote 
by V^ the union of those V, for which = 3, and by V^ the union 
of those V, for which = 4. Now (, the number of 3's in V, equals 

c = c (3) + c (4) 
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where (^ 3 \ respectively denotes the number of 3's in V^ 3 \ respec- 
tively . Therefore 

Prob[C > a\V\] < Prob[C (3) > ^\V\] + Prob[C {4) > ^\V\] 

Z z 

We will separately estimate the two terms in the right-hand side. The 
last one is the easiest. By definition of all spin values equal to 3 in 
are contained in the islands listed above. Each island contains < 8 
squares and has an energy > 4. So the energy per square is > 1/2. If 
£( 4 ) > | V | a/2, then the energy of the configuration is > a|V|/4. Note 
that E 3 y > 1 because the configuration = 3 has zero energy. So the 
probability of the configuration rj is less than or equal exp(— f3a\V\/4). 
As the total number of configurations does not exceed 4' y ', 

^(C (4) >f|V|)<4^exp(-^|V|) (4.5) 

which, when af3 > 41og4, tends to zero as \V\ — > oo. 

The rest of the proof is an estimate of ^yiC^ > a 1^1/2) • 
Note that C (3) < \V (3) \, and hence 

/4,,(C (3) >%\V\)< 4,(1^1 > f \V\) = A. (4.6) 

where Z a is the partition function of those configurations i] e 7Zv,3 for 
which \V^\ > § \V\. 
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It is easy to see that 



Z <Ve-" |r| " (12) - (12) f4 7l 

r 

where the sum runs over those polycontours T, for which |\^ 3 )| > ^ \V\ 
(r = 0, for which V {3) = V is included). 

Let us now estimate "B^y. Consider only those configurations 77 G 
7Zy,3 which also belong to 7Zy, 4 , where V' is V without squares touching 
9V. Then we have 

■3 > _-3/3|9V|-(12) o\ 

From the cluster expansion (or by direct counting), 

E m > e -W\ E M (4.9) 

for some constant h > 0. Therefore, 



-3,V 



> e-^lsfj (4.10) 



for some constant d (depending on (3). So we obtain 



-/3|r|^(i2) -(12) 
d\dv\ ^yp>"4,yw 



Za ^ d\av\ ^ r e 



where the sum over T is as in (J4.7j) . In the right-hand side use that for 
every £x G 7^ (3) 4 and for every £ 2 £ ^y(4) 45 there exists a ( 6 7?^ 4 
that coincides with £1 on V^ 3 -* and with £2 on . Hence, 



= (12) ^ -(12) w (12) (A , \ 

"4,y — "4,y< 3 ) "4,y< 4 ) 
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and thus 

"(12) 

I^<e d ^J2^ m ^g^ (4.13) 
" 3,y r ' =, 4,y( 3 ) 

Now use the inequality (|4.4|) in order to obtain the following estimate 

(4.6) <(4.13) < e «fl^l^ e -/'|r| e /|avW| e -c|v(3)|e-^ (414) 

r 

But evidently |aV^^ 3 > | < \dV\ + \T\, and W 3 ) > a|F|/2, so we have the 
estimate 

eid+mv \ e -i*c\v\e-wj2e-™ (4.15) 

r 

where /3' = (3 — f . To estimate the last sum is standard combinatorics: 

r x&X \ 7 3a;,|7|>12 J 

where X is the set of all vertices — ends of edges inside V, \X\ < \V\. 
A classical Peierls' estimate gives 

<T m < H'e- 1313 ' < He^ (4.17) 

~f3x,\-y\>12 

for some constants H = H' exp[13/]. Since 

1 + He~ 13t3 < e He ' 13 " 

we obtain the estimate 

(4.6) < e (d+f)\9V\ e -±ac\V\e-W e H\V\e-W (41g) 
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If (3 > \og[2H/ac], this product tends to zero as \ V\ — > oo, |dV|/|V| — > 
0. That together with (4.5) concludes the proof. 

Remark 4.6. To prove that the density of l's or 2's goes to zero is even 
easier: Let N be the number of l's and 2's in V, N = iV (3 - 4) + N [T] 
where is the number of l's and 2's in V (3) U V {A) and is the 

number of l's and 2's in [r]. Observe that < 2|T|, and the l's and 
2's in U are contained in islands. For the islands, the energy 
per square is > 1/2. The energy of a configuration is equal to |T| plus 
the energy of the islands. Hence, if N > cx\V\, then that energy is not 
less than a|V|/2. We can then do as in (4.5) to prove that 

l4,v( N ^ a \ V \) as|V|^+oo 
whenever a (3 > 2 log 4. 

The same proof as above remains valid for the ensemble TZy,4- In fact, 
it is valid modulo straightforward modifications, for general boundary 
conditions as well. The proof for the ensemble with general bound- 
ary conditions £ only needs small modifications concerning the islands 
touching the boundary. 
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